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ABSTRACT 

Until recently, the methods of geodetic positioning on the earth were limited to a precision of 
roughly one part in 10 6 . At this level of precision, the earth can be regarded as a rigid body since 
the largest departure of the earth from rigidity is manifested in the strains of the earth tides which 
are of the order of one part in 1 0 7 . Long baseline interferometry is expected to routinely provide 
global positioning to a precision of one part in 1 0 8 or better. At this level of precision, all parts of 
the earth’s surface must be regarded as being, at least potentially, in continual motion relative to the 
geocenter as a result of a variety of geophysical effects. 

This paper discusses the general implications of this phenomenon for the theory of the earth’s rota- 
tion and focuses particular attention on the question of the measurement of the “earth’s rotation 
vector” on a deformable earth. 


♦Present! address until September 1, 1979: Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, 
Pasadena, California, 91103. 
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INTRODUCTION 

Almost 1 0 years ago, a paper entitled “Geodetic Ties Between Continents by Means of Radio Tele- 
scopes” [Jones, 1969] appeared in the journal Canadian Surveyor . The author, Harold Jones of the 
Geodetic Survey of Canada, reported on what was probably the first geodetic application of the 
technique of long baseline interferometry (LBI). Jones’ solution for the 2400 km baseline vector 
between antennas at Prince Albert, Saskatchewan, and Algonquin Park, Ontario, yielded a “correc- 
tion” of 30m + 20m to the value previously determined by ground-based survey methods. 

Baseline vectors of comparable magnitude are now routinely determined by LBI to an accuracy of 
±10 cm to ±100 cm. This represents an improvement of two orders of magnitude in 10 years and 
stands as a testimonial to the determination and skill of those who have worked to develop the tech- 
nique over the past decade. So dramatic has been this development that LBI measurements of the 
near future, with expected accuracies of ±1 cm to ±5 cm between continents, coupled with other 
known geophysical effects, such as continental drift, threaten to render obsolete our fundamental 
definitions of the body-fixed and space-fixed coordinate frames which serve as a basis for defining 
positions on the earth and their motion relative to inertial space. In discussing this issue, it is useful 
to present a brief historical description of the evolution of these concepts, for the present situation 
is best understood in terms of its past history. This past history is in turn often useful in indicating 
the direction of the next logical step in the development. 

The rotation of the earth or the earth’s orientation in space is specified by the time dependent 
transformation between a set of body-fixed basis vectors (e 2 , , £ 3 ) conceived of as being 

“attached” to the earth and a set of space-fixed basis vectors (E 1 , E 2 , E 3 ) conceived of as being 
“attached” to space. In general, the time dependent transformation can be represented as 

e ; (t) = T y (t) El (1) 

where the transformation matrix T^ft) expresses all that is known about the earth’s rotation. (Ex- 
planatory Supplement to the Astronomical Ephemeris and the American Ephemeris and Nautical 
Almanac, 1960; Mueller, 1969). Our understanding of T^ft), which has progressively improved 
historically, can be usefully classified into six distinct eras. 


Ancient Era 

Since ancient times, mankind has identified two great circles on the celestial sphere. The first was 
the equator whose location on the celestial sphere was defined by the diurnal rotation of the earth. 
The second was the ecliptic whose location on the celestial sphere was defined by the annual mo- 
tion of the brightest object in the sky - the sun. 

According to the ancients, these two great circles were fixed on the sky, providing reference direc- 
tions and a celestial coordinate frame with which to develop elaborate sky maps populated with 
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mythical characters. The zodiac, which is referenced to the ecliptic and the equator, dates from 
these early times. 

If these early conceptions were to be updated into modem language, they could be described by the 
introduction of two sets of basis vectors with a common origin at the center of the earth consisting 
of: 


(a) a space-fixed set of basis vectors, Ej, E 2 , E 3 , such that E t is contained by the equator and 
the ecliptic at the ascending node, E 3 is oriented toward the pole of the ecliptic, and E 2 
makes up a right-handed orthogonal triad; and 

(b) a body-fixed set of basis vectors, ij , e 2 , e 3 , such that e 3 is parallel with the earth’s rota- 
tion axis, ej is contained by the intersection of the equator and the prime meridian, and 
e 2 makes up a right-handed orthogonal triad. 

With the exception of the sun, the moon, and the five planets observable with the naked eye which 
were deified for their persistent anomalous behavior, the ancients described all motion by spinning 
the earth at a uniform rate about an axis fixed in the earth and fixed in space. The mathematical 
transformation between the space-fixed basis vectors and the body-fixed basis vectors could be 
represented by the introduction of a time dependent spin matrix Sjj(t). 

e.(t) = Syft) Ej (2) 

The quantities dSjj/dt were considered to be constants since the earth’s rotation was held to be 
uniform. 


The Classical and Medieval Era (500 B.C. - 1600 A.D.) 

During this period, the precision with which astronomical observations could be made steadily im- 
proved, as did the mathematical models required to make astronomical predictions of eclipses, oc- 
cultations, and other periodic celestial phenomena. During this era, no truly dynamic theories of 
celestial motions were introduced, and the algorithms on which calculations depended were, for the 
most part, highly geometrical. The principal advance in the knowledge of the earth’s rotation which 
occurred in this era was the discovery by Hipparchus, circa 150 B.C., of the precession of the 
equinoxes. 

Following Hipparchus’ discovery, previously held conceptions of the space-fixed basis vectors Ej 
were retained; however, it was necessary to redefine the space-fixed basis vectors as being the E ; 
referenced to some celestial pole and equinox of a particular epoch. The body-fixed basis vectors 
e. retained their previous definitions. 
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In this era, the practice of realizing the body-fixed basis vectors by assigning numerical coordinates 
to observatories fixed on the earth was developed. A similar practice of realizing the space-fixed 
basis vectors by assigning numerical coordinates to stars at a particular epoch led to the develop- 
ment of extensive star catalogues. 

With the discovery of precession, the mathematical transformation between the space-fixed basis 
vectors and the body-fixed basis vectors required the introduction of the time dependent precession 
matrix P..(t). 

S,(t) = P u (t)S Jk (t)E k (3) 


The quantities dS^/dt were still held to be constants. 


The Newtonian Era (1600 A.D. - 1800 A.D.) 

The Newtonian era was characterized by the rapid development of Newton’s dynamical theory and 
its extensive application to geodynamics and problems of the earth’s rotation. Around the year 
1680, Newton offered an explanation of the precession of the equinoxes and predicted the oblate 
figure for the earth on the basis of his theory of gravitation. At about the same time, Newton and 
Huyghens published a predicted value for the polar flattening of the earth. 

These important geophysical developments were being matched by corresponding dramatic astro- 
nomical discoveries. In 1718, the astronomer Edmund Halley announced the discovery of stellar 
proper motion, and the researches of the astronomer James Bradley, circa 1728, led to his discovery 
of both stellar aberration and astronomical nutation. This was followed in 1735 by the geodetic 
surveys in Finland and Peru under Pierre Bouguer of the French Academy which verified Newton’s 
prediction of polar flattening and confirmed the theory of precession and nutation. 

The discoveries of the Newtonian era placed the definition of the space-fixed coordinate frame on a 
firm dynamical footing. As a result of the discovery of the periodic phenomenon of nutation, it be- 
came necessary to reference the space-fixed basis vectors Ej to the mean equator and ecliptic of a 
particular epoch, averaging out the small periodic effects of nutation on the location of the equator. 

In addition, the discoveries of the Newtonian era made it necessary to realize the space-fixed basis 
vectors E, by assigning coordinates, not to the stars themselves which were no longer regarded as 
being fixed, but to locations on the celestial sphere known as “mean stellar positions of epoch” 
which were related to the actual stellar positions by: 

(a) the removal of the effects of proper motion since the epoch, and 

(b) the averaging out of the periodic effect of annual aberration. The relatively constant con- 
tribution to the aberration arising from the effects of the eccentricity of the earth’s orbit 
were retained in the stellar mean position. 
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The body-fixed basis vectors retained their classical and medieval definition and continued to be 
realized by assigning coordinates to observatories fixed on the earth. The meridian through Green- 
wich was widely adopted as the prime geographic meridian at this time. 

With the discovery of nutation, it became necessary to introduce a third time dependent matrix, the 
nutation matrix N^ft), into the mathematical transformation between the space-fixed basis vectors 
Ej and the body-fixed basis vectors e r 

e i (t) = P ij (t)N jk (t)S k£ (t)E s . (4) 

The quantities dS^/dt were still held to be constant. 


The Euler/Kelvin Era (1800 A.D. - 1900 A.D.) 

The Euler/Kelvin era was characterized by two major developments in the theory of the earth’s rota- 
tion. The first was inspired by the work of Leonard Euler who in 1765 published his famous trea- 
tise on the dynamics of rigidly rotating bodies. Euler showed that for an oblate spheroid such as 
the earth, the rotation axis, the angular momentum axis, and the figure axis will not coincide in gen- 
eral and furthermore, that in the event that they do not coincide the rotation axis cannot remain 
fixed within the body. 

As a consequence of Euler’s work, a search was begun circa 1840 for evidence of the effects of 
Eulerian motion on the earth’s rotation. This search persisted for 50 years before achieving success, 
and at one time or another involved the active participation of Peters, Bessel, Maxwell, Kelvin, New- 
comb, Kiinster, and Chandler (Munk and MacDonald, 1960). 

The second major development of this era was inspired by the work of Kelvin (W. Thompson, 1863) 
who in 1863 was the first to introduce into theoretical geodynamics considerations of the deforma- 
bility of the earth. Kelvin’s seminal work on tidal deformations of the earth stimulated the growth 
of an entire branch of geodynamical research concerned with the general geodynamical effects of 
earth deformations which included the work of Darwin, Lamb, Love, and others and which is being 
vigorously pursued to this day. 

These two developments were wed in a remarkable manner when, in 1891, the work of Kiinster and 
Chandler finally revealed the elusive Eulerian motion of the earth with a period 40 percent larger 
than that predicted by Euler’s rigid body theory. This gross discrepancy between the theoretical 
and observed period of the wobble was reconciled the following year by Newcomb, who showed 
how the yielding of the deformable earth to the centrifugal forces of the wobble reduced its effec- 
tive dynamical ellipticity and lengthened the wobble period to the observed value. 

The Euler/Kelvin era ended in 1900 with the establishment of the five observatories of the Inter- 
national Latitude Service (ILS). 
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The developments of the Euler/Kelvin era produced no essential modifications of the definition of 
the space-fixed coordinate frame spanned by the basis vectors E r The equator of rotation and the 
equator of figure were now regarded as distinct and moving relative to each other. However, their 
relative motion was such that their mean positions respectively coincided with the equator of angu- 
lar momentum, and since the basis vectors E ; were already defined in terms of the mean equator 
and ecliptic of epoch, the space-fixed basis vectors Ej retained to their previous Newtonian defini- 
tions and realizations. 

The discoveries of the Euler/Kelvin era profoundly altered the definitions of the body-fixed coordi- 
nate frame spanned by the basis vectors e r In the first place, since the rotation axis could no longer 
be regarded as a body fixed axis, it could no longer serve as a reference axis for the e 3 basis vector. 
Consequently, e 3 was redefined as being parallel to the earth’s figure axis. However, as a conse- 
quence of the demonstration of the earth’s tidal deformability, the figure axis was not fixed within 
the earth. It was necessary to finally define the e 3 body-fixed basis vector as lying in the direction 
of the mean figure axis of the earth which is related to the instantaneous figure axis by averaging 
out the effects of deformations due to tides and polar motion. This averaging was actually carried 
out by ILS observatories over the years 1900 to 1905. 

The e x body-fixed basis vector was then defined as lying in the planes of the mean equator of figure 
and the meridian through Greenwich which has been adopted internationally as the prime meridian. 
The body-fixed basis vector e 2 was defined to complete a right-handed orthogonal triad. 

The body-fixed basis vectors e 3 were now realized in practice by assigning coordinates, not to the 
observatories themselves which were no longer regarded as being fixed, but to their mean locations 
on the earth which were related to their actual locations by: 

(a) averaging out latitude and longitude variations due to polar motion, and 

(b) averaging out tidal displacements. 

With the discovery of Eulerian motion, it became necessary to introduce a fourth time dependent 
matrix, the wobble matrix W^(t), into the mathematical transformation between space-fixed basis 
vectors E 4 and body-fixed basis vectors e^ 

g i (t) = P..(t)N jk (t)S ke (t)W Em (t)E m (5) 

In this stage of the development of the theory of the earth’s motion, the instantaneous rotation axis 
is taken as an observable reference axis and the orientation of the earth (spanned by Sj) in space 
(spanned by E } ) is expressed by (see figure 1): 

(a) First, orienting the reference axis in space, accomplished by the matrices Pjj(t), Ny(t); and 

(b) Second, orienting the earth relative to the reference axis, accomplished by the matrices 

s«(t). W ij( t). 
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The matrix S^Ct), which describes the orientation in azimuth of the earth about the reference in- 
stantaneous rotation axis, no longer has the property that the quantities dS^/dt are constants. This 
is merely a consequence of the kinematical choice of reference axis and the fact that the chosen 
reference axis is itself moving within the body of the earth as is explained by Woolard (1953, pp.27, 
163-165). However, while they are no longer constants, the kinematical quantities dS^/dt are pre- 
cisely known from the dynamics, and consequently, the elements of the matrix S;.(t) are precisely 
predictable. 

The matrix WL(t) which describes the location in the body-fixed frame of the instantaneous axis of 
rotation can be decomposed into two components: 

(a) a predictable contribution W7(t) due to the effects of the gravitational torques of the sun 
and moon and often referred to as the “dynamical variations of latitude,” and 

(b) an unpredictable part W? (t) due to a variety of internal geophysical effects. The unpre- 
dictable part of Wjj(t) is determined observationally and published in arrears. 
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Early Twentieth Century (1900 A.D. - 1960 A.D.} 

The era extending from roughly 1900 A.D. to 1960 A.D. was characterized by a rapid and unprece- 
dented expansion of geophysical research and international cooperation. Many important advances 
in instrumentation and observation techniques occurred which had bearing on the measurement of 
the earth’s rotation. Some of these include the development of the zenith telescope (1912), the 
PZT (1915), the quartz crystal clock (1950), the dual rate moon camera (1954), and the atomic 
clock (1958). 

Technological development resulted in a rapid sequence of startling geophysical discoveries in all 
branches of geophysics, oceans, atmosphere, and “solid” earth including the earth-moon system and 
“near space” as well as the two principal fields associated with the earth - the gravitational field and 
the electromagnetic field. The principal impact of this era on the theory of the earth’s rotation 
came from the discovery by Stoyko in 1936 that the earth’s rotation described by the matrix S^(t) 
contained a variable component which was not predicted by the existing theory. This effect dis- 
covered by Stoyko was largely due to the seasonally varying atmospheric winds. 

Throughout this era beginning with the work of the ILS and its successors, the IPMS and the BIH, 
the body -fixed basis vectors Bj were referenced: 

(a) e 3 to the mean figure axis of the earth as established by an average over the interval 
1900 to 1905, and 

(b) Sj to the meridian of the mean observatory and the plane of the mean figure axis of 
1900 to 1905. 

The “mean observatory” is a fictitious observatory occupying the meridian to which the time scales 
UT1 and UT2, published in the forms of UT1-UTC and UT2-UTC, are referred. In general, the 
meridian of the mean observatory does not contain the transit circle of the observatory at Green- 
wich but usually passes to the east or west of Greenwich by a few milliseconds of time depending 
on the particular mix of observatories and their respective weightings which have contributed to the 
data base from which the time scales UT1 and UT2 have been deduced. The concept of the “mean 
observatory” has been introduced to allow for variations in the least squares adjustment procedures 
by which UT1 and UT2 are obtained (Mueller, 1969). 

During this era, the body-fixed basis vectors B ; were realized by assigning coordinates to mean ob- 
servatory positions which have been obtained from the actual observatory positions by: 

(a) averaging out latitude and longitude variations due to polar motion, 

(b) averaging out tidal displacements, and 

(c) allowing for the effects of the least squares adjustment procedure. 
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Throughout this era, the space-fixed basis vectors £L retained their previous Euler/Kelvin era defini- 
tions and realization. 

The mathematical transformation from the space-fixed basis vectors E. to the body-fixed basis vec- 
tors Sj is still expressed as it was in the Euler/Kelvin era, 

( 6 ) 

except that now both S^t) and Wy(t) consist of a combination of predictable and unpredictable 
contributions. In the case of S-ft), the predictable contribution constitutes the majority of it with 
the unpredictable portion appearing as a small perturbation. In the case of WU(t), the unpredictable 
contribution constitutes the majority of it with the predictable portion appearing as a small 
perturbation. 


The Present Day Era (1960 - ) 

The year 1 960 is identified here as the epoch marking the beginning of the present-day era of the 
study of earth rotations, for in the few years surrounding 1960, two scientific discoveries occurred 
which at the time seemed entirely unrelated to each other and to the question of earth rotation but 
may ultimately prove to be of fundamental importance in advancing our future understanding of 
earth rotation. 

First was the geophysical discovery of what has come to be known as “global plate tectonics.” This 
“discovery” actually spanned about a decade of time and, like most discoveries which involve a large 
number of participants, was the product of, and to some extent anticipated by, earlier work. Never- 
theless, this decade of discovery spanned the years 1955 to 1965 and can be regarded as being made 
up of two somewhat distinct intervals: 

(a) 1955 to 1960 during which the palaeo-magnetic evidence of large scale polar wander also 
confirmed the hypothesis of continental drift, and 

(b) 1960 to 1965 during which the conception of “sea floor spreading” with rigid tectonic 
plates in relative motion at their boundaries emerged and was confirmed by the studies of 
the polarity reversals of the geomagnetic field and the magnetic striping on the ocean 
floor at the mid-ocean ridges. 

These discoveries destroyed the validity of the concept of the “mean location” of an observatory 
and along with it destroyed the rigorous validity of the operational procedure by which the body- 
fixed basis vectors Bj are realized. According to the theory of plate tectonics, the observatories are 
moving secularly relative to each other and relative to the bulk of the earth in quasi-random direc- 
tions and at rates which vary between 1 and 1 0 cm/year. Under these circumstances, the concept of 
an observatory’s mean position has no inherent meaning, for it is a function of the averaging inter- 
val and epoch chosen. 
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The discovery of plate tectonics injected a certain lack of rigor into the observational procedures 
which had developed historically for measuring earth rotation and created a vague discomfort in 
the minds of geophysicists who worried about such matters. However, it had little impact on the 
day to day operations of earth rotation measurement, because the effects being manifested by con- 
tinental drift and plate tectonics were below the detectability threshold of the instruments of the 
instruments of the 1 960’s. This situation changed dramatically with the development of the long 
baseline interferometer. 

The second event was the astrophysical discovery in 1 960 of quasi-stellar objects often referred to 
as QSO’s or quasars. By 1963, the quasars were known from their red shifts to be among the most 
distant objects in the universe. The long baseline interferometer was developed in 1967 specifically 
to study the radio emission from these objects. 

It was the long baseline radio interferometer which linked the separate discoveries of sea floor 
spreading and quasars to the question of the earth’s rotation. This was due to the fact that, in the 
operation of a long baseline interferometer, there was for the first time the potential to make meas- 
urements of the earth’s rotation so accurately that the effects of the secular deformation rates im- 
plied by the theory of global plate tectonics could be detected. Confronted with this fact and the 
continued development of the long baseline interferometer as a device for making regular earth 
rotation measurements, it is rapidly becoming necessary to develop rigorous operational procedures 
and definitions which explicitly incorporate into the measurement of each rotation the secular de- 
formability and continuing fracture of the earth’s crustal structure. 


Strain in Deformable Bodies: The Separation of Rotations and Deformations 

We begin by considering two mass elements of a continuous deformable body which before defor- 
mation occupy points Pj and P 2 given by position vectors and"? 2 respectively and separated by a 
vector dx. The geometry of this arrangement is illustrated in figure 2. The body is then subjected 
to a deformation which carries the mass elements at Pj and P 2 through displacements u, and u 2 to 
points Qj and Q 2 respectively. If the deformation field u(r) is a continuous and differentiable func- 
tion of position, which is equivalent to asserting that no fracture has occurred, then for points P : 
and P 2 separated by an infinitesimal vector dx, we have 

u 2 = Uj + Vu * dx! (7) 

The relative displacement vector Au = u 2 - Uj is given by 

Au = Vu . dx (8) 

where 


( 9 ) 
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The tensor Au can be split up into its symmetric and antisymmetric parts 


where 


1 1 „ 
Vu = — e.. + — J2.. 
2 IJ 2 « 


( 10 ) 


is symmetric and where 


e :; = 


3uj 3uj 


1J 3Xj 3Xj 


3Uj 3Uj 
3Xj 3Xj 


( 11 ) 


( 12 ) 


is antisymmetric, which allows the components of relative displacement Am to be written 


AUj = J e y d Xj + j S2 y d Xj . 


(13) 
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Like all antisymmetric tensors, £2^ possesses a dual tensor given by 


e ijk ^jk 


where is the alternating tensor defined by 

r. 


e Uk = < 


[+1 ijk cyclic 123 
0 ijk not all distinct 
^T ijk noncyclic 123. 

It can easily be shown from the property of that 

n ij = e kij n k’ 

and, as a result, the components of relative displacement can be written 

Au i = 7 e ij dx j + ! e kij n k dx r 
By introducing the axial vector 3 given by 

ej = e.. k n jk e., 

we see that equation (16) can be written in coordinate free notation as 


— ► 1 ^ 1 

Au = — e . dx £2 X dx. 

2 2 


( 14 ) 


(15) 


(16) 


(17) 


(18) 


We see from this result that, in general, the relative displacement of two mass elements of the me- 
dium separated by an infinitesimal vector dx can always be decomposed into the sum of a relative 
displacement Au D , 


Au r 


e . dx, 


(19) 


arising solely from the effects of deformation (should any be present) and a relative displacement 
Au„ , 


Au r = - -AX dx, 


( 20 ) 
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arising solely from the effects of rigid body rotation (should any be present). In this way, the ef- 
fects of deformational and rotational contributions to the relative displacement field are separable. 

The strain tensor e whose elements are e^ has the property that if ds 2 is the squared distance be- 
tween the mass elements at the points P 1 and P 2 before deformation and ds' 2 is the squared dis- 
tance between the same two mass elements after deformation, then 


5(ds 2 ) = ds' 2 - ds 2 = e^ dXj dXj. (21) 

The rotation tensor £2 whose elements are £2^ has the property that the axial vector £2 obtained 
from its dual gives the amount of rigid body rotation suffered by the separation vector dx as a re- 
sult of the displacements. 


An alternate and, in this instance, quite useful approach to the problem [Brillouin 1964, pp. 287 
ff . ] is to adopt the point of view that the mass elements of the medium are assigned unchanging co- 
ordinates whose numerical values are preserved under the material transformations which accom- 
pany the displacements. This requires the coordinate surfaces to be associated with the material 
medium and to be carried along and deformed with it as it moves. These deformed coordinate sys- 
tems generally will be curvilinear. 

In the usual treatment of this problem, the mass element at position 7 and coordinates x p related by 

7= x iei + x 2 e 2 + x 3 e 3 , (22) 

gets displaced to a position 7' = 7+ u(r) with coordinates X ; such that 


and 


X i = X i + Uj 


(23) 


r ’ Xl e 1+ X 2 e 2 + X 3 e 3 


(24). 


or equivalently 


7' = (Xj + u^ej + (x 2 + u 2 )e 2 + (x 3 + u 3 )e 3 . (25) 

In the alternate approach, the mass element at position 7 and coordinates x ; gets displaced to a posi- 
tion?' =7+ u(r) with coordinates Xi, such that 


and 


X i = *i 


r ' = X, e' + X„e; +X,e 


l l 


2 2 


3 3 


(26) 

(27) 


or equivalently 
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where the basis vectors e' , e' 2 , locally spanning the X coordinate frame are necessarily changed 
from the basis vectors e 3 , e 2 , e 3 locally spanning the X coordinate frame in order to accommodate 
the fact that "r# ?/ 

Without loss of generality, we can begin by assuming that the basis vectors e 3 , e 2 , e 3 span a local 
Cartesian coordinate frame. The relationship between the coordinates Xj referred to the basis vec- 
tors e t , e 2 , e 3 , and the coordinates X ; referred to the basis vectors e' , e 2 , e 3 can be established by 
considering two mass elements at position r } and ? 2 separated by the infinitesimal vector dx, where 

dx = r 2 -r r (29) 

As a result of the^ displacements and the motion of the mass elements, the vector dx is transformed 
into the vector dX: 


dX = ? 2 - rj = ? 2 + u(7 2 ) ~ ?i - u (?j ) (30) 

or 

dX = dx + du (31) 

where 

du = u(? 2 ) - u(?j). (32) 

Now, in component form, equation (31) gives 

dXj = dXj + dUj, (33) 

and the chain rule of differentiation gives 

dX. dx. 9u. dx- 

— - = — i . (34) 

dx k dx k d Xj dx k 

Recognizing that dXj/dXj = 5^ where 5 y is the Kronecker delta, and that by the definition of equa- 
tion (26) dX k = dx k , we obtain the relationship that 


dXj = 


6., + 

ik 


du. 


9x b 


dX k . 


The use of equation (10) in equation (35) gives 


dXj = 


1 1 „ 
5., + — e.., "1 T2., 

ik 2 lk 2 lk 


dX,. 


(35) 


L 


(36) 
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which states that the X and X coordinate systems are locally transformed, one into the other, by 
the combination of (a) a rigid body rotation £2 jk and (b) a deformation e jk . 

If we proceed to compute the squared distance ds 2 between the two mass elements after the dis- 
placements on both the X coordinate frame and the X coordinate frame, these two expressions must 
be equal since ds 2 is a scalar quantity whose value is independent of the choice of coordinates. In 
the X coordinate system, which has been assumed without loss of generality to be a local Cartesian 
coordinate frame with metric tensor g^ = 5^, then 

ds 2 = 5.; dX, dX.. (37) 

A 1J 1 J 

The equivalent expression in the X coordinate frame is obtained from equation (36) as 

ds | = 5 ij 5 ik + \ e ik + \ n ik 5 je + \ e je + \ n i2 d *k A \ ( 38 > 

which, using the antisymmetric property of reduces finally to 

ds|=(5 ij +e ij )dX i dX j (39) 

or 

ds| =g tj dXj dXj. (40) 

A comparison of equations (37) and (39) shows that a consequence of the theoretical development 
on which the displaced mass elements are allowed to retain their coordinate values is that the strain 
tensor e^ is simply absorbed into the geometric tensor g^ of the deformed X coordinate system. 
Whereas the basis vectors e x , e 2 , e 3 are all orthogonal and of unit length, the basis vectors e j , e ' 2 , e 3 
are nonorthogonal and no longer of unit length. This is illustrated in figure 3. 



Figure 3 . 
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The metric tensor gy of the X coordinate frame is given by definition as 

|e'| 2 le'Me^l cos e[e^ |e'| |e^| cose[e^ 

le^le'l cos e^e^ le^l 2 |e^| |e^| cose' 2 e' 3 

|e'||ej| cose'e' |e^| \e' 2 \ cose^ le' | 2 

So, we see from equations (39), (40), and (41) that 


Sij = 


(41) 


|8j| - \/l + e n 

I e 2 1 = ^ + e 22 


le' | = Vl +e 


33 


(42) 


and that 


cos = 


'12 


V 1 + j V 1 + e 


22 


cos e^e^ =■ 


'13 


\/ 1 + e j j V 1 + e 


33 


(43) 


cos e^ = 


'23 


Vl + e 22 V 1 + e 


33 


If the | e y | « 1 then these relationships reduce to 


1 


le;i= 1 +— e 


l S 2l- 1 + T 6 22 


2 11 

1_ 

2 


- 1 + 7 e 33 


( 44 ) 
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and to 


e i e 2 


7 r 

— - e. , 
2 12 


e i e 3 


e 2 e 3 


IT 

2 13 
TT 

= T " 623 


(45) 


The entire preceding analysis is based on the components of the separation vector dx and the com- 
ponents of the relative displacement vector Au. The direction of the separation vector dx is arbi- 
trary and a reversal of the direction of dx simply reverses the sign of Au. The quantities in this 
theory have analogs in long baseline interferometry. The analog of the separation vector dx is the 
interferometer baseline vector (f separating the two antennas. The analog of the relative displace- 
ment vector Au is A the change in the baseline vector which is observed to occur between two suc- 
cessive measurements of j?. 


Next we shall attempt to apply these procedures to the analysis of long baseline interferometry data 
and to question of the measurement of earth rotation on a deformable earth. While these proce- 
dures are potentially useful in this context, they will also have certain limitations. 


First, the above geometrical procedures will not give rise to tensor quantities when applied to the 
whole earth. Tensor quantities can only be defined in the neighborhood of a point or, for pairs of 
points, separated by infinitesimal separation vectors dx. In long baseline interferometry, we will be 
concerned with relative displacement between points separated by thousands of kilometers. The 
analysis procedure described above will produce quantities which have many of the properties of 
tensors but, because they refer to pairs of points separated by a finite distance, will not possess the 
general transformation properties of tensors. In particular, these quantities will possess the mathe- 
matical properties of tensors only in Cartesian coordinate frames. To distinguish them from true 
tensors, we will refer to them as tensor analogs. 


Second, the deformation field of the earth is not continuous and differentiable everywhere. Frac- 
ture is occurring on and within the earth in a complex, largely unknown pattern. This fact does not 
affect the validity of these procedures as applied to the earth insofar as we restrict ourselves to the 
discipline of kinematics or geometry. However, entering the realm of dynamics or geophysics and 
attempting to relate strain tensor analogs to stress tensors will be a risky procedure. Where there is 
fracture, there can be local strain without local stress. Without considerably more geologic knowl- 
edge, the complexity of the pattern of fracture affecting baselines of thousands of kilometers in 
length would make the inference of stress from strain suspect. The two possible exceptions to this 
might be: 
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(a) The case of strain tensor analogs measured on very short baselines from a few tens to a 
few hundred kilometers in length, in which case, the material supporting the interferom- 
eter can either be considered homogeneous and free of fracture or can be mapped in suf- 
ficient geological detail to adequately account for the effects of fracture; and 

(b) The case of a mean strain tensor analog obtained by a global average of many local strain 
tensor analogs, in which case, the random effects of geologic faulting affecting the local 
measurements have been averaged out and the residual mean strain tensor analog de- 
scribes that global strain field which is deforming the earth coherently as though it were 
the neighborhood of a point. 


The Measurement of Earth Rotation on a Deformable Earth 
Basic Definitions and Procedures 

We presuppose a global network of interferometer baseline vectors /3 k k = 1, 2, 3, ... N forming a 
polyhedron as shown in figure 4 which will be used for the purposes of: 

(a) defining the body -fixed coordinate frame, and 

(b) measuring earth rotation or equivalently defining the transformation from the body-fixed 
coordinate frame to the space-fixed coordinate frame. 


63 



Figure 4 . 
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These definitions must be operable in the face of the fact that as a result of the discoveries of plate 
tectonics, each point of the earth’s surface must be regarded as possessing its own rotation axis. 
While points on the same tectonic plate might be expected to have nearly identical rotation axes as 
a consequence of Euler’s theorem on rigid body motion on the surface of a sphere [McKenzie, 
1972] , different rotation axes will be associated with points on different tectonic plates. Further- 
more, when the relative vertical motions, which occur within the plates and which can be as large as 
1 cm/year, are also considered, it is clear that in general the concept of a rotation axis must be rigor- 
ously regarded as only applying locally on the earth. 

Each of the time dependent interferometer baseline vectors /^(t), k = 1, 2, 3, ... N will possess 
space-fixed components B k (t), i = 1, 2, 3, 

£ k (t) = B k (t) Ej + B k (t) E 2 + B k (t) E 3 , (46) 

and body-fixed components b k (t), i = 1, 2, 3, 

$ k (t) = b k (t) e t +b k (t)e 2 + b k (t)e 3 . (47) 

The transformation from space-fixed components to body-fixed components expresses all our 
knowledge of earth rotation and is currently expressed by equation (6) as 

f(t) = Wj(t) SjTft) Nj,(t) pj m (t) B* (48) 

where the superscript “k” serves merely to identify the baseline and does not obey the Einstein 
summation and range convention. The superscript “T” denotes a matrix transpose. 

Repeated measurements of the baselines p k (t) at times . . . t - t , t„ t .. t yield rela- 

r r v 7 m— 2 m-1 m m + 1 m + 2 J 

tive displacement vectors A k (t m ) given by 

2' t ( , „) = 3 k (t m )-3 k (49) 

The relative displacement vectors A k (t m ) will possess space-fixed components D k (t), i = 1, 2, 3, 

2 l (t m )=D!|(t m )E, +D5(t n ,)E 2 +D5(t nl )E 3 (50) 

and body-fixed components d k (t m ), i = 1, 2, 3, 

^ k ( t m) = d l( t m) S l +d 2< t m) g 2 +d 3( t m>®3 ( 51 > 

also related one to the other, according to our present knowledge of earth rotation, by equation (6) 

> - Wj(t m ) S,T t (t„ ) Nf,(t m ) Pj m (t m ) Dj, (t. ). 


(52) 
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The operation of the long baseline interferometer and the subsequent reduction of the data finally 
provide us with the quantities b k (t m ), dK (t m ) i = 1 , 2, 3, k = 1, 2, 3 . . . N which together with the 
epochs t m , m = 1 , 2, 3 . . . constitute the basic data with which we must achieve our two stated 
objectives. 


We can begin creating analogs to continuum mechanics by defining, for the k th baseline, a dimen- 
sionless quantity cj^ given by 



d k b k 


b k b k 

u s 


(53) 


which will serve as the analog of the tensor Vu of equation (9). The tensor analog cK is expanded in 
its symmetric and antisymmetric parts. 


c k = — e k ; + — fi k (54) 

u 2 IJ 2 1J 

where 

e k - = [cf. + cK] = — ^ — [d k b k + d k bh (55) 

J J J b k b k 

is the strain tensor analog for the k th interferometer baseline and where 

ft k = = [c k - c k ] = — 1- [d k b k - d k b k ] (56) 

1J J1 i 1 J J 1 

b 5 b S 

is the rotation tensor analog for the k th interferometer baseline. It can readily be shown that these 
definitions permit us to write 


d f = T e » b f + T^ b f- <57) 

Equations (54) through (57) are the analogs of equations (10) through (13) for the case of contin- 
uum mechanics. 

To illustrate this procedure, we consider the case of an interferometer baseline p 1 * with body-fixed 
components b j' , i = 1, 2, 3, given by x, y, z and a relative displacement vector A k with body-fixed 
components d k , i = 1 , 2, 3, given by Sx, 5y, 5z, then: 
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c Er 


x 2 + y 2 + z 2 


eh = • 


x 2 + y 2 + z 2 




X 2 +y 2 + z 2 


x5x y5x z6x 

x5y y5y z5y 
xSz y5z z5z 

2x6x y5x+x5y z5x+x5z 
x5y+y5x 2y5y z5y+y5z 
x5z+z5x y5z+z5y 2z5z 

0 y5x-x5y z5x-xSz 
x5y-y§x 0 z5y-y§z 
x5z-z5x y5z-z5y 0 


(58) 


(59) 


(60) 


To show that eh possesses the properties of a strain tensor insofar as the k th baseline is concerned, 
we compute the variation in the squared distance s k2 between the two observatories defining the 
k th baseline which occurs as a result of the relative displacement. Before relative displacement 


and after relative displacement 


The difference 5(s k2 ) is given by 


s k2 = bh b k , 

s k2 = (b k + d k ) (b k + d k ). 


(61) 

(62) 


5(s k2 ) = (b k + d k ) (b k + d k ) - b k b k 


or 


5(s k2 ) = 2d k bh 


1 +0 


2dh b k 


(63) 


neglecting terms of the order of d/b « 1. Equation (63) can be written 

5(s k2 ) = 


dhb k + d k b k 




bh b k = eh b k b k 


(64) 


which is the analog of equation (21) for the case of continuum mechanics. 

A further demonstration that eh possesses the properties of the strain tensor for the k th baseline is 
provided by considering the special case of an earth subjected to a uniform dilatation. In this case, 
the observatories at position ? 1 and ? 2 which make up the k th interferometer baseline p k suffer dis- 
placements 5? 1 and S? 2 given by 
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Sfj = Afj 


Sr 2 = Ar 2 


where 3 A is the magnitude of the dilatation. The relative displacement vector A k is given by 

^ k =A(? 2 -r 1 ) 


and will have body-fixed components d k 


d k = Ab k 

where b k are the body-fixed components of the k th baseline vector j3 k = ? 2 The strain tensor 
analog for the k th baseline in this case is 

e 5 = ^ [b 'VW <«> 


and the trace of eK, T r (e k ), in this case is 


T r (e k ) = 2A. (66) 

This result should be compared to the well known result from continuum mechanics that 

T r ( e ij) = 20 (67) 


where Q is the dilatation. 

To show that £2 k . possesses the properties of a rotation tensor insofar as the k th baseline is con- 
cerned, we compute the rigid body rotation c5 which is imparted to an interferometer baseline vec- 
tor p as a result of a relative displacement A. The geometry of this situation is illustrated in figure 5 
with the body-fixed basis vectors shown for reference. The rotation vector co can be expressed as 

(3 = ( 68 ) 

where £2 is the magnitude of the angle between 0 and |3 + A. Neglecting terms of the order of A/j3 
— ► — ► 

where j A| « lp|, we can express the area of the shaded triangle in two equivalent forms and equate 
them to obtain 


Y£2|^| 2 = -j |/? x A|. 


(69) 
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From this it follows that 


a. I?x 4l . 

\$\ 2 

Since 

0X A 

"" $X A| ’ 

equations (70) and (71) give 

- $x£ 

CO = ’ 

i?i 2 



In the body-fixed frame, 





I ii mi ■ ■■ ■■ ■■ in am ■■■■ 
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and so 


= e jn h ^ 

' i> s b t 


The antisymmetric tensor dual to the axial vector cu is 


“ij “ e k ij “k 


or 


OJ.. = 
IJ 


e kij e k8m b S d m 


Using the properties of the alternating tensor, this can be written 


( 73 ) 


(74) 


or 


b 0 d_ 

2 m 

w.. = (5-„ 5. -5. 5 0 ) 

i] v iS jm lm ^ 


(75) 


“ij 


b i b j - b j d i 
b n b n 


(76) 


A comparison of this result with the definition in equation (56) of the rotation tensor analog £2ji 
for the k th baseline shows that 


< 77 > 

As a consequence of this, we see that the body-fixed components £2* of the vector of the k th 
interferometer baseline occurring as a result of the rotation relative displacements A k are obtained 
from the rotation tensor analog S2 k j by 


=- 


€■ 


im n 


f2 k 


(78) 


The Body-Fixed Coordinate Frame 

From observations on a global network of N interferometer baselines, we can define 



k = l 


N 

£w k e k 

k= 1 


( 79 ) 
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and 



k = 1 


(80) 


where W k k = 1, 2, 3, . . . N are appropriate weighting factors. Ey is the weighted mean global strain 
tensor analog and Ay is the weighted mean global rotation tensor analog. The strain tensor analog 
e k - and rotation tensor analog £2K on each interferometer baseline can then be expressed as the sum 
of the above global means and a “local” residual e k and coK respectively specific to the k th baseline. 


e k = E.. + e k 
u ij u 


\ + ^5 


(81) 

(82) 


where, by definition, the weighted global means of the local residuals vanish. 


1 


N 



k= 1 


L 

v = 1 


W k e k = 0 


(83) 


co k = 0. (84) 


k= 1 

If the global distribution of interferometer baselines is arranged to representatively sample the 
earth’s crust, then Ey might be expected to reveal the deformation properties of the earth as a 
whole. Such measurements might be expected to reveal a number of global processes which are 
expected to be producing coherent deformation of the whole earth. 

Earth expansion or contraction would be revealed by a significant departure of T. (Ey ) from zero. 
Overall earth expansion or contraction is associated with a number of outstanding problems in geo- 
physics such as: 

(a) the problem of the internal evolution of the earth and chemical phase changes possibly 
associated with the growth or decay of the solid inner or liquid outer cores, 

(b) the problem of the effect on the earth’s figure of the secular spin-down and the question 
of the finite strength of the mantle and the persistence or otherwise of a “fossil” equa- 
torial bulge, and 

• 

(c) the problem of possible cosmological effects due to a non-zero value of G/G causing secu- 
lar changes in the earth’s size. 
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The off-diagonal elements of Ey would be associated with global shear deformation and would re- 
veal tidal-like effects on the earth. The eigenvectors of Ey would indicate the directions of princi- 
pal global stress as evidenced by crustal observations. Their evolution with time might reveal infor- 
mation about processes in the earth’s interior. Of particular interest would be a comparison of the 
evolution in time of the eigenvectors Ey as measured by long baseline interferometry and the evolu- 
tion in time of the principal axes of inertia of the earth as evidenced by its rotational dynamics and 
the orbits of satellites. 

For certain geophysical applications such as the study of plate tectonics, it may prove useful to 
introduce a material body-fixed coordinate frame x, y, z spanned by basis vectors e,, e 2 , e 3 de- 
fined in the sense that for coherent global deformations Ey, the body-fixed coordinate surfaces, 
deform along with the global motion of the material. This is easily accomplished using the observed 
mean global strain tensor analog Ey according to the procedure described above in the section on 
continuum mechanics. 


At some epoch t Q , the earth’s internal state is adopted as an arbitrary standard with Ey = E?. = 0 
and the material body-fixed basis vectors e 2 , e 2 , e 3 coinciding with the basis vectors e 3 , e 2 , e 3 de- 
fined in equation (6). The material coordinates x, y, z are therefore identical to the Cartesian co- 
ordinates x, y, z at the epoch. For times t > t later than the epoch, the material body-fixed basis 
vectors e. , e 9 , e, deform along with the mean global motion and are given by equations (44) and 
(45) as 


' 5 ll= l+ 7 E ll 


le 2 l=M-- E 22 


(85) 


and 



e ie 2 =tt/ 2-E 12 
e 2 e 3 - tt/2 - E 23 


( 86 ) 


The coordinate system x, y , z spanned by basis vectors e 3 , e 2 , e 3 might offer certain advantages in 
studies of plate tectonics as the coordinate velocity, or changes in the coordinates of points on the 
earth referenced to this coordinate frame, will be due only to “local” tectonic effects. The effects 
of coherent global deformation are removed automatically by deformation of the coordinate frame. 
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The strain tensor analog Ey represents mean strain in the global network of interferometer base- 
lines shown in figure 4. The residuals given in equation (81) represent local deformations of the 
k th baseline relative to the global network and are due to the effects of regional tectonic motion. 
The eK will contain the information about regional tectonics on the earth including both inter- and 
intra-plate geologic processes. 

If the global distribution of interferometer baselines is arranged to representatively sample the 
earth’s crust, then the mean global rotation tensor analog Ay and its dual axial vector A with body- 
fixed components Aj given by 


A i = e ijk A ij 


(87) 


represent a net mean rotation of the network of interferometer baselines shown in figure 4 relative 
to the body-fixed basis vectors Sj , e 2 , e 3 of equation (6). If the global network of interferometers 
is to be used to define the body-fixed coordinate frame spanned by 8j , e 2 , e 3 , it is a contradiction 
for such a coordinate frame to be rotating relative to itself. Hence, the body-fixed basis vectors e,, 
e 2 , e 3 must necessarily be defined in such a way that Ay vanishes. Since Ay can be regarded as an 
infinitesimal rotation residual to the transformation of equation (6), the most direct way of defin- 
ing 8j , e 2 , e 3 so that Ay vanishes is to include the effects of the rotation of Ay directly into the 
transformation of equation (6) from the space-fixed to the body-fixed frame. This can be done by 
rewriting the transformation as 


®i(t> = [*u - \(t)] P jk (t> N ke (t) (t) W m n (t) E n (88) 

However, this is probably not the most useful procedure to follow since the measured values of 
Ay(t) contain, in addition to the unknown effects of global deformability, the effects of errors in 
the models Py, Ny, etc. for the known effects of precession, nutation, spin, and wobble. By observ- 
ing the time dependence of the measured values of Ay, it will be possible to decide between the 
alternatives of modifications to the known transformation matrices Py, Ny, Sy, Wy and the addition 
of yet another transformation matrix Dy representing the rotational effects of global deformability. 

Historical precedent indicates that the discovery of new geophysical or astronomical phenomena has 
been accompanied by the addition of a new matrix to the transformation from the space-fixed 
frame to the body-fixed frame. However, since the models for precession, nutation, etc. are also 
expected to contain errors at the level of measurement anticipated by future long baseline interfer- 
ometry, we should be prepared to exercise both options. 


The residuals coji given in equation (82) represent local rotations of the k th baseline relative to the 
global network and, like the residuals Cy of equation (81), are due to the effects of regional tec- 
tonic motion. The coK, like the eK , will contain information about regional tectonics on the earth 
including both inter- and intra-plate geologic processes. 


136 


RADIO INTERFEROMETRY 


The Space-Fixed Coordinate Frame 

The present set of space-fixed basis vectors E 1 , E 2 , E 3 are dynamically defined by the orbit of the 
earth about the sun and, while they play an essential role in the formulation of the theory of the 
orbital dynamics of the earth, they play no essential role in the formulation of the theory of the ro- 
tational dynamics of the earth. The basis vectors Ej , E 2 , E 3 have served both purposes in the past 
as a matter of convenience. However, since translational (orbital) and rotational motions of the 
earth are dynamically independent, no great difficulty should arise if the two motions were referred 
to different space-fixed coordinate frames - should it prove convenient to do so. It is apparent that 
the application of long baseline interferometry to problems of earth rotation would be greatly facil- 
itated by referencing the space-fixed basis vectors directly to the cosmic radio frame of the quasars. 
Besides being the most remote objects in the observable universe, hence expected to possess very 
little proper motion, the choice of the quasars to define the space-fixed frame offers some impor- 
tant measurement advantages. 

The delay observable r k is a scalar quantity related to the inner product of a unit vector ? in the di- 
rection of the radio source and the baseline vector (P 1 

r k = -s.fr. 
c 

The delay observable can be written in two equivalent forms 

= — s 1 
c 1 

r k = — c. |3 ik 
c 1 


(89) 

(90) 

(91) 


where s 1 , s i5 are the contravariant and covariant components of the unit source vector s and j3 ,k , j3 k 
are the contravariant and covariant components of the k th baseline vector j? k both referred to a set 
of, as yet unspecified, basis vectors the usual choice for which is E x , E 2 , E 3 or ^ , e 2 , e 3 . 

If we introduce a space-fixed coordinate frame whose fundamental directions are defined by three 
selected quasar radio sources s(l), s(2), s(3), then it is clear that the baseline vector /3 k has two 
representations given by 


fr = (3 k s 1 + 0 k s 2 + j3 k s 3 (92) 

and 

fr =0 lk Sj +|3 2k s 2 +J3 3k s 3 (93) 

where s 1 , s 2 , s 3 are the contravariant quantities or so-called basis 1 -forms [Misner, Throne, Wheeler, 
1973, pp. 53 ff. ] and Sj , s 2 , ? 3 are the covariant quantities or so-called basis vectors spanning the 
space-fixed frame appropriate to radio sources s(l), s(2), s(3). 
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An analysis of the operation of a long baseline interferometer [Thomas 1972, Cannon 1978] shows 
that the delay observable yields the covariant components of the baseline vector as shown schemat- 
ically in figure 6. The stopped fringes of the incoming wave fronts from the radio source provide a 
physical realization of the basis 1 -forms which define the contravariant quantities s 1 ,? 2 , s 3 . In the 
space-fixed coordinate frame defined by radio sources s(l), s(2), s(3), the contravariant basis vectors 
or 1 -forms s 1 , s 2 , s 3 have components 


s 1 = [ 1 0 0 ] 

s 2 = [ 0 1 0 ] (94) 

s 3 = [ 0 0 1 ] 

by definition. From equation (90), the baseline vector has representation 

j f k = C7^ S 1 + C7^ S 2 + CT^ s 3 (95) 

where the covariant space-fixed components of $ k given by 

= cr l 

j3* = cr^ (96) 

are thus measured directly without reference to source coordinates. 



Figure 6. 
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It is traditional in geodesy to measure contravariant components of baseline vectors - the infinitesi- 
mal displacement vector being defined as the archetypical contravariant vector. To obtain the con- 
travariant space-fixed components {3 lk of the baseline vector /3 k from the covariant space-fixed com- 
ponents j3 k given by the delay observable and equations (96) requires knowledge, of the metric 
tensor gy of the nonorthogonal space-fixed coordinate frame of the quasars s(l), s(2), s(3) shown 
schematically in figure 7. The contravariant components (3 lk are given by 


where 
and where 


P ik = g ij Pf 


g^gy 1 


Sij = 


1 coss(l)s(2) coss(l)s(3) 

coss(2)s(l) 1 cos s(2)s(3) 

coss(3)s(l) cos s(3)s(2) 1 


(97) 

(98) 


(99) 


This procedure does not require knowledge of absolute source coordinates but merely the relative 
angles or the direction cosines between the quasars s(l), s(2), s(3). This measurement can also be 
made very accurately by long baseline interferometry. 



We now consider the case where observations are made of a radio source s(J), J = 4, 5, 6 . . . which 
is not one of the three sources s(l), s(2), s(3) defining the fundamental directions in the space-fixed 
frame. The delay observable on the source s(J) and baseline (3 k is r k (J) where 

r k (J) = — s(J) • p k (100) 

c 

and where s(J) is a unit vector in the direction of the radio source s(J). s(J) can be written in terms 
of the funamental contravariant basis vectors? 1 , s 2 , s 3 as 
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s(J) = s 1 (J)s 1 +s 2 (J)s 2 +s 3 (J)s 3 (101) 

where 

Sj (J) = cos s(l)s(J) 

s 2 ( J) = cos s(2)s(J) (102) 

s 3 (J) = cos s(3)s(J) 

Equations (102) indicate that any of the quasar radio sources s(J), J = 4, 5, 6 . . . can be observed to 
obtain the space-fixed covariant components of the interferometer baseline providing the direction 
cosines of the unit vector s(J) relative to the three chosen fu ndam enta l spac e-fixed directions given 
by sources s(l), s(2), s(3) are known. The angles s(l) s(J), s(2) s(J), s(3) s(J) can be measured accu- 
rately by long baseline interferometers. This suggests that a cataloguing procedure in terms of the 
direction cosines to three chosen fundamental sources instead of right ascension and declination 
might be usefully adopted in the near future. 

It is possible that the three chosen radio sources s(l), s(2), s(3) will exhibit proper motion relative 
to the cosmic space-fixed radio frame defined by the ensemble of quasars. Should this be the case 
it will present no great difficulty to the above procedure. The ensemble of quasars would then be 
chosen to define the space-fixed frame and to possess a Maxwellian velocity distribution and to ex- 
hibit no mean rotation. With this assumption, it would remain only to tie the positions of the 
chosen sources s(l), s(2), s(3) to the ensemble of quasars by periodic measurements of the angles 
s(l) s(J) sC2)~s(J) st3)~s(J), J = 4, 5, 6 . . . N where N could be of the order of 100 to correctly ac- 
count for any proper motion they might exhibit. 

Finally, it may perhaps be useful in the light of the general theory of relativity to retain and to dis- 
tinguish between two sets of space-fixed basis vectors. The first set? 1 ,? 2 ,? 3 described above span 
what could be called a cosmic inertial frame. The second set E t , E 2 , E 3 to be defined dynamically 
by a relativistic ephemeris of the solar system span what could be called a local inertial frame. The 
theory of general relativity allows for these two inertial frames to be relatively rotating [Weinberg, 
1972, pp. 239 ff]. 
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APPENDIX 

The Definition of the Mean Global Strain Tensor E- 
and the Mean Global Rotation Tensor A- 

The definition of the mean global strain tensor Ey and the mean global rotation tensor Ay as given 
in equations (79) and (80) 



N L 

k=l 


N 

I w 

k=l 


k k 
'ij 




k=l 



k=l 



can be shown to satisfy two desirable objective criteria. These are discussed below. 


E- - the Mean Global Strain Tensor 

The definition of E ;j . minimizes, in a 2 weighted least squares sense the net departure between the 
squared baseline length variations 5s which would be predicted assuming a mean global strain 
E.j and the individual squared baseline length variations 5s k2 observed by the long baseline 
interferometer. 


r k k , k , k 
os = e . . b . b . 

iJ i J 


<5s k = E. . b k b k 
^ 1 J 


— ^2 2 

Net weighted squared departure between predicted 5s and observed 5s k squared baseline length 
variations is D 2 where 
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Iw k (6s k -^i k2 ) 2 

k=l 
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D 2 = Iw k [(e k .b k b k ) 2 - 2e k . E b%% k b k + (E b^ k ) 2 ] 
k-1 1 ij i J lj Jim i j £ m £m £ m 


A choice of E 2m which minimizes D 2 is given by 


3D 

3E_ 


= 0 


pq 


This gives a condition 


N 


l » k l-2<^ -5T 2 b M b K + 2(E ij b ^ b j ) W L(E *>* b 
=1 t pq pq 

J iW k (- 2ek . bjbjbjb k 6 £p 6 mq + ZE-.btt^ 6 £p 6 mq 

+ 2 v^) = 0 


;) = o 


= 0 


f W k f-2e k . + 2E..]b^bJb k b k = 0 

^ [ ij ijj 1 j p q 


Since this condition should be fulfilled independently of the particular choice of baseline net work 
we conclude that we shall in general require 



0 
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or 



1 N 

N k l w 

I w k=l 
k=l 


k k 

ij 


Ajj - the Mean Global Rotation Tensor 

The definition of A.j minimizes, in a weighted least squares sense, the magnitude of the net depar- 
tures between the mean global rotation vector A ; and the rotation vectors A? 1 of the individual base- 
lines observed by the long baseline interferometer. 


A. 
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E i£m 
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G i£m 



The magnitude of the net weighted squared departure between the mean global rotation vector A ; 
and the observed rotation vectors £2^ of the individual baselines is e 2 

N 

6 2 = y w k (n k -A.) (fi k -A.) 

, i i i i 

k=l 


but since 


ft k -A . = e.„ (fi k -A. ) 

l l i£m £m £m 


we have 


® 2 = E ipq tn L~ A £m 1 [£i pq- A pq 3 

K. 


From the properties of the alternating tensor this can be written 
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e 2 = l w fc 6 0 6 [nj - a, ][n k -A ] 
£p mq pq pq 
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e 2 = l w k [n k -a. ][n k -A ] 

^2. £m £m Jim £m 


6 2 = l w k [fl k 2 -2A ft k +A 2 ] 
, S £m £m £m Jim 
k=l 


A choice of A„ which minimizes 6 2 is given by setting 


This leads to the condition 


« , on. . on. 

y w k f-2 — — o k +2A — = 0 

L . 1 M W £ni £m 3A J U 
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y W [-2fif 6. 6 +2A. 6 n <5 1 = 0 

_ 1 £m £p mq £m £p mq J 


f w k [-2ft k + 2A ] = 0 
k=l pq pq 


The solution to this equation is simply 
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